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A study on the behavior of off-site AC power failure recovery times at three nuclear plant sites is 
presented. It is shown, that power law is appropriate for the representation of failure frequency- 
duration correlation function of off-site power failure events, based on simple assumptions about 
component failure and repair rates. It is also found that the annual maxima of power failure 
duration follow Frechet distribution, which is a type II asymptotic distribution, strengthening our 
assumption of power law for the parent distribution. The extreme value distributions obtained are 
used for extrapolation beyond the region of observation. 
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1. Introduction 

Estimation of off-site power failure characteristics is important for the safe design and operation of nuclear power 
plants. The emergency power supply requirements are based on the ability to forecast the maximum credible loss of 
off-site power (LOSP) duration during the life time of the plant. The observed data is available only over a period 
of less than two decades and scarce on long duration failures. We address two questions: first, what is the nature of 
the parent distribution, i.e., the form of frequency-duration correlation function and second how to extrapolate for 
CZ2 ■ failure durations longer than the observed maximum. 

First we show, based on simple assumptions about the off-site power supply system component failure and repair 



. rates, that the frequency-duration correlation function is a power law of the form, y cx x . Next, it is established that 



the distribution of annual maxima, of observed failures follow Type II asymptotic distribution (Frechet distribution, 
see Johnson 1970). This result strengthens the premise, that the parent distribution for power failure duration- 
correlation is power law in character. Extreme value analysis is also used to make extrapolations for more number of 
years than observed, for instance to estimate the most probable loss of off-site power duration in say, 50 or 100 years. 



> ■ 

2. Off-site power failure frequency-duration function 

The off-site AC power at a plant may be lost due to supply failure from the power grid or because of plant centered 
equipment failure like station transformer, feeders, breakers etc. The grid failure could be a minor grid disturbance 

C*~) ■ or grid collapse. The grid power may also fail due to severe environmental conditions. Baranowsky, (1988) has used 
Weibull functions of the form (Aje^ - " 4 *' ') for representing the frequency of off-site power failure events of each type, 
exceeding a given duration. Although the fits are good, it would be better if some physical basis could be established 
for the fitted function. Further considering the scarcity of data, instead of a series of stretched exponentials a single 
function would be preferable. It is postulated here that there is a mixture of power failure rates Xi with a certain 

•i-H | probability density g(\i) as A^ G {Xi\i = L.iV}. Similarly for the repair rates ^ G — 1..N}. That is, there are 
different kinds of equipment each with a specific Xi and /x^. Then the probability of observing a failure of duration 
r-| ' t' > t is 

. £ : P(t' >t)=J2 9(k) exp(- M it)<2 (1) 

><■ * 

where e^ Alit - ) is the probability that the failed component is not restored in time t, assuming constant repair rate 
/ij. Since more frequent failures are generally repairable in shorter times compared to rare failures which take more 
time to set right, it is assumed that fi cx A, with this assumption in equation^ and going over to the continuous limit 
we get, 

P(t'>t)= f g{\)e { - cXt) d\ (2) 
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Where c > is a constant of proportionality. The form of g is not known, however it is known (Johnson, 1970) that 
if gamma distribution 

g(X) = ae(- aX \aX)( k -V /T(k) 

for constants a, k > 0, is assumed for g, then Pareto distribution is obtained as the cumulative distribution function 
F(t) = 1 - P(t' > t). That is the solution of equation H is 

For large t, P(t) goes like P{t) - (9/t) k ; k > 0; 9 > 0; t > 0, where 6 = a/c and 

F(t) « 1 - (0/t) k (4) 

This family includes Cauchy, t and F distributions. The function P(t) is plotted for data observed at two plant sites, 
over a period of 15-20 years, in Fig. 1 and 2. The plots show P(t) versus LOSP time t in Log-Log scale. The linear 
fits are shown alongside, which confirm the power law character of the distribution. In Fig. 1 there is considerable 
deviation from straight line for the first few points. This is due to the fact that for small values of t, i.e., for —t << 1, 
P(t) is virtually constant, and the power law dependence manifests fully for large t. The LOSP frequency-duration 
correlation function has the same time dependence as it will be different from P(t) only by a constant factor. 

3. Extreme value analysis of LOSP duration 

The distribution of maximum of m observations of the random variable t, distributed as F(t) is, 

lim F(t) m = G{t) 

711 >GO 

The non-trivial limit distribution G(t) is obtained by appropriate scaling of the variable t 7 

lim F((t - b m )/9 m ) m = G(t) 

m ►oc 

when F(t) is of the form of equation^! the asymptotic limit, with b 7n = is 

G(t) = e-WV" (5) 

which is one of the three possible asymptotic limits, that is, a type II extreme value distribution (Fisher et al., 1928 
and Gumbel, 1954) also known as Frechet distribution. For finite m equation [S] is written as G m (t) = e~( 6 ' m / t ) fc . 

Extreme value analysis of LOSP data collected in (Theivarajan 1999, Marimuthu et al., 2000 and Kumar et al., 
2002), is done as follows. The annual maximum ij for each location is arranged in ascending order as, t\ < t2 < ••• < £jv- 
The plotting points F(i) are calculated as 

F(i) = (£-0.3)/(iV + 0.4) 

(See Karl Bury, 1999), where N is the number of years over which the data is collected. The F(i) approximate the 
median of the distribution free estimate of the cumulative distribution function. When — ln(— ln(F(i))) are plotted 
against ln(i), where t is the power failure duration, straight lines are obtained as depicted in Figs. 3, 4 and 5, as 
expected from equation [S] The exponents k obtained from the slope, is shown in the respective figures. 
The distribution of maximum in any n years will be 

G n (t)=G m {t) n =e-lW (6) 

where 

n = 9 m n Wk) 



The most probable maximum in n years of observation, obtained from equation [S] is, 

f = 0(_2*_)(V*) 
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The most probable maximum t is plotted in Fig. 6, for the three locations. The relatively higher value of most 
probable maximum for TAPS and FBTR cases is due to higher incidence of recorded on-site power distribution 
equipment failure. From Fig. 6 we can infer that, for instance a 20 h LOSP extrema is likely in TAPS site in 15 y, 
whereas events of such duration are expected to occur only in 50 years and 100 years respectively for the FBTR and 
MAPS sites. These conclusions are based on the data collected from these sites and results are indicative only. 

4. Conclusion 

The nature of the relationship of off-site power failure duration and its frequency is studied and it is found to 
have power law dependence. Plausible physical basis for the observed behavior is given. The power law nature of 
asymptotic behavior is confirmed by performing an extreme value analysis of the data. The extreme value distribution 
is also used to extrapolate beyond observed LOSP durations. 

It is interesting to note that, power laws have been observed in a variety of natural and man made settings like, 
frequency distribution of words (Zipf, 1949), distribution of incomes, earthquake magnitudes and recently in internet 
topology (Faloutsos et al., 1999), to name a few. Based on the intimate connection between power laws and the ideas 
of self organized criticality (Bak et al., 1987), it can be said that off-site power distribution system seems to evolve into 
a critical state, where failures of longer duration are not infrequent, as one would expect for exponential distribution. 
Although it is convenient not to consider cutoffs to failure duration from a theoretical perspective, it would be better 
if methods could be devised to include cutoffs considering the finitcness of the system. 
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FIG. 1: P(t) versus loss of off-site power duration t for site Kalpakkam-MAPS 
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FIG. 2: P(t) versus loss of off-site power duration t for site Kalpakkam-FBTR 
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FIG. 3: Log Log of Fi versus Log of maximum annual loss of off-site power duration for site Kalpakkam-MAPS 
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FIG. 4: Log Log of Fi versus Log of maximum annual loss of off-site power duration for site Kalpakkam-FBTR 
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FIG. 5: Log Log of Fi versus Log of maximum annual loss of off-site power duration for site Tarapur-TAPS 
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FIG. 6: Most probable maximum LOSP duration versus time of observation in years for the three nuclear plant locations 



